A=[ 1  2  1 ; 1  0  -1 ; 3  8  3];
b=[0 0 -2]';
x=A\b
A*x